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Abstract 

We study the Lagrangian dynamics of systems of point vortices and passive particles 
in a two-dimensional, doubly periodic domain. The probability distribution function of 
vortex velocity, p^i has a slow- velocity Gaussian component and a significant high- velocity 
tail caused by close vortex pairs. In the limit for N ^ oo, pi^j tends to a Gaussian. However, 
the form of the single-vortex velocity causes very slow convergence with A^; for A^ ^ 10^ 
the non-Gaussian high-velocity tails still play a significant role. At finite A^, the Gaussian 
component is well modeled by an Ornstein-Uhlenbeck (OU) stochastic process with variance 
(Jat = ^A^ In N/2-K. Considering in detail the case A^ = 100, we show that at short times 
the velocity autocorrelation is dominated by the Gaussian component and displays an 
exponential decay with a short Lagrangian decorrelation time. The close pairs have a long 
correlation time and cause nonergodicity over at least the time of the integration. Due to 
close vortex dipoles the absolute dispersion differs significantly from the OU prediction, 
and shows evidence of long-time anomalous dispersion. We discuss the mathematical 
form of a new stochastic model for the Lagrangian dynamics, consisting of an OU model 
combined with long-lived close same-sign vortices engaged in rapid rotation and long-lived 
close dipoles engaged in ballistic motion. From a dynamical-systems perspective this work 
indicates that systems of dimension O(IOO) can have behavior which is a combination of 
both low-dimensional behavior, i.e. close pairs, and extremely high- dimensional behavior 
described by traditional stochastic processes. 
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1. Introduction 

Advection of passive and active tracers is a crucial component in many geophysical 
processes: ozone transport in the stratosphere, pollutant dispersal in the atmosphere and 
the ocean, plankton and salinity transport in the ocean. Additionally, understanding the 
relationships between Eulerian and Lagrangian statistics is necessary for interpreting data 
provided by ocean floats and atmospheric balloons. 

Standard Lagrangian approaches to transport in fluids are based on the use of 
either stochastic models or chaotic advection in low-dimensional dynamical systems. 
Lagrangian stochastic models are generally based on assuming that transport is dominated 
by unstructured Gaussian random velocity fluctuations, and they are most successful when 
the system under study has an extremely high-dimensional phase space. The standard 
example is Brownian motion which describes the irregular movement of microscopic 
particles in systems with dimensionality of the order of Avogadro's number, 0(10^^). In 
this case, the details of the deterministic description are irrelevant. 

Chaotic advection, on the other hand, is based on a fully deterministic description of the 
phase-space dynamics. It has traditionally been applied to systems with just a few excited 
degrees of freedom, with very simple spatial structures, and with periodic or quasi-periodic 
temporal dynamics.!' B i Thus, the traditional application of stochastic models is to 
highly turbulent flows without strong coherent structures, while low-dimensional chaotic 
advection is appropriate for flows dominated by a few large-scale waves, vortices, or other 
structures.!' i 

Geophysical turbulence, however, does not fully belong to either of the above categories. 
One of the main difficulties encountered in studying Lagrangian advection in geophysical 
turbulence is, in fact, the complex space-time structure of the flow. The joint effects of 



- 3 - 



rotation and stratification often induce the presence of energetic coherent structures that 
contain the majority of the enstrophy of the system. ME Intense jets such as the Gulf 
Stream act as partial barriers to transport, while coherent vortices such as ocean mesoscale 
eddies and the stratospheric polar vortex can trap particles for long times. 
Previous studies of barotropic turbulence have shown that coherent vortices play an 
important role in the advection of Lagrangian tracers, inducing characteristic signatures 
that cannot be captured by simple stochastic models.0' 

In this work, we further explore the properties of advection in flows dominated by 
strong coherent vortices, and, in particular, we consider an ensemble of many point vortices 
in two spatial dimensions. Systems of points vortices capture many of the features of 
two-dimensional turbulent flows.S EE 0' El In geophysically relevant situations the 
number of energetic coherent structures is neither as large as Avogadro's number, nor as 
small as the number of degrees of freedom in simple models of low-dimensional chaotic 
advection. Systems of point vortices with ~ 0(10™), m small, are thus used here as 
another simplified paradigm: they have an intermediate phase-space dimensionality, and 
are appropriate for describing advection in vortex dominated flows. As we show below, 
the dynamics of systems of many point vortices displays properties of both stochastic 
models and low-dimensional chaotic advection, and provides a useful bridge between 
high-dimensional stochastic models and low-dimensional chaotic advection. 

The remainder of this paper is organized as follows. Section 2 contains an introduction 
to the dynamics of point vortices. Section 3 addresses the central limit theorem for 
point-vortex systems and considers the statistical properties of ensembles of point vortices. 
In Section 4 we discuss the statistical properties of long time-integrations of systems of 
100 point vortices, and study the behavior of both the vortices and of passively advected 
particles. In section 5 we study single-particle dispersion in point-vortex systems and 
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compare with the dispersion of an ensemble of randomly moving particles described by the 
Ornstein-Uhlenbeck (OU) stochastic process.! Finally, in Section 6 we present conclusions 
and perspectives, and discuss a possible alternative stochastic model for describing advection 
in point-vortex systems. 



2. Point- Vortex Dynamics 

Point vortices are singular solutions of the Euler equations in two spatial dimensions. 
The dynamics of an ensemble of point vortices is described self-consistently by a system of 
equations that take the non-canonical Hamiltonian form 

dt di/i ' dt dxi ' 

where Xj = (xj, is the position of the i-th vortex with constant circulation Fj, H is the 
Hamiltonian 

^({x.}) = - E ^G(x„x,), (2) 

is the number of vortices, and the form of the Green function G depends on the boundary 
conditions.ii The positions Xi and yi play the role of non-canonically conjugate variables, 
the number of degrees of freedom is equal to the number of vortices, and the dimensionality 
of the phase space is twice the number of vortices. The number of independent conserved 
quantities of (p and (0) also depends on the boundary conditions, but it is always finite 
and small. Hence, an assembly of more than a few vortices behaves chaotically. 

Passively advected particles are easily incorporated in point-vortex systems. A point 
vortex with zero circulation, Fj = 0, is advected by the other vortices but has no influence 
on the velocity of any other vortex. Such passively advected particles will be referred to here 
as passives, and particles with Fj 7^ as vortices. Recently, the motion of passive particles 
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in point- vortex systems has been the subject of several investigations.^ H These studies 
have shown that point vortices are surrounded by finite-size islands of regular Lagrangian 
motion, where passives can be trapped for very long times. This implies the existence 
of long-time non-ergodicity in passive particle motion, where averages over ensembles of 
different advected particles differ from time averages over single-particle trajectories. In 
turn, this may lead to long lasting differences between Eulerian and Lagrangian averages 
and the possibility of non-Brownian (anomalous) dispersion. In the following, we explore in 
detail some of these issues. 

In this work we use periodic boundary conditions because it is the only domain that has 
the needed properties. A closed domain has a finite maximum separation and is thus not 
suited to study long-time dispersion properties. This leaves the infinite plane, the periodic 
plane (2-torus), and the sphere. The infinite plane is unsuitable because motion is not 
homogeneous with a finite number of vortices. On a sphere, there is no unique way to count 
the number of times a particle travels around the domain, and thus the sphere does not 
allow a satisfactory definition of particle dispersion at long times. Thus the only boundary 
condition which meets our requirements is that of the periodic domain. In addition, most 
simulations of homogeneous turbulence use these boundary conditions. 

The Green function G for point vortices on a periodic domain with length 27r may be 
written aS 

^ ^ [ cosh {xj ~Xj- 27rm) - cos {y^ - yj) \ (x, - xjf 
G(x„x,j- 2. inj^ cosh i2nm) ) 2n ' 



r?i=— oo 



The function G can be shown to be periodic in x and y, and invariant under the 
transformation x ^ y. The velocity of the i-th vortex resulting from (|T]), (j^), and (|^) is a 
sum over the velocities induced by each of the other vortices: 
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Ui = yT- y -sm{yi-yj) ^ 

^ ^t^^ cosh(xi - Xj - 27rm) - cos(?/i - yj) ' 

y.^yY- y sm(x^ - Xj) ^ 

i=i ^ m=-oo cosh(?/, - ?/j - 27rm) - cos(xi - Xj) ' 

where u = {ui,Vi) = {dxi/dt,dyi/dt). The dynamics on the periodic domain has two 
invariants, the components of the hnear momentum, corresponding to translation parallel 
to the X and y axes; = Y.f=i ^iVi and Py = Y^iLi ^i^i- Angular momentum is not an 
invariant because the periodic boundary conditions break the rotational symmetry. The 
invariants Px and Py are independent only if the total circulation J2 is zero.@ In this 
case, the motion is chaotic for > 3. 



3. Point-Vortex Configurations 

In this section we discuss the properties of instantaneous configurations of vortices. 
We restrict ourselves to configurations with zero total circulation and an equal number of 
positive and negative vortices. 

First consider the velocity induced at a point xq = {xo,yo) by a single vortex at xi. 
In the limit r = |xi — xqI — the velocity (Q) asymptotically approaches the velocity 
generated by a single vortex on the infinite domain, 2r/r.Q At larger distances from the 
vortex, the domain periodicity modifies this simple form of the velocity. In practice, the 



*Note that with the choice (0) for the Green function, the timescale differs from that 
usually used on the infinite plane by a factor of An; hence the velocity from a close vortex is 
here a factor An greater than the expression usually used on the infinite domain.^! 
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velocity on the periodic domain departs significantly from tlie infinite plane result only for 
separations greater than approximately one sixth of the domain. 

Next, consider the probability density function (pdf) for the velocity resulting from a 
single vortex, pi{u). This has been studied in detail by Min et al in the case of the infinite 
domain.ii If the position of the vortex is chosen at random from a uniform distribution, 
then the high velocity limit of pi is easily calculated. Since the probability of having 
separation between r and r + dr is 27rrdr, and since for small r the velocity scales as 1/r 
then for large |u|, pi{u.) ~ j9(r(|u|)) |(ir/(i|u|| ~ l/|up. This scaling for large velocities is 
independent of the boundary conditions, as the boundaries are irrelevant at sufficiently 
small separations. The detailed form of the distribution for smaller |u| does, however, 
depend on the specific form of u for a periodic domain, and it depends on the direction of 
u as well as its magnitude. 

Since we are interested in the dynamics of systems of several point vortices, we now 
consider the velocity induced at a point by a system of vortices with |rj| = 1. This 
velocity is merely the sum of the velocities due to each vortex separately. If the vortices 
are randomly placed, then the velocity is the sum of random numbers, and its pdf 
]9Ar(u) is determined by pi(u). For functions pi which decay sufficiently rapidly, the central 
limit theorem applies and pn approaches a Gaussian. Most statements of the central limit 
theorem require that the variance of pi be finite. Here, since pi ~ l/l'^P ^^e variance 
/ |up]9i(|u|)(i|u| diverges logarithmically. However, this divergence is sufficiently slow that 
the central limit theorem still applies and a.s N ^ (yo, pn does become Gaussian. For more 
details, see the Appendix, Min et al^l, and Ibragimov and Linnik.S 

For large but not infinite A^, the velocity pdf of an ensemble of vortices displays 
interesting properties. For distributions with finite variance, the convergence to a Gaussian 
is quite rapid. For example, if pi were constant over a finite region and zero elsewhere. 



- 8 - 



then is quite close to Gaussian even for = 5. However, for distributions with 
slowly diverging variance, such as the case here, the convergence to a Gaussian is slow 
and is called "non-normal" B For non-normal convergence, the variance of the sum scales 
differently than the normal case. In the Appendix we show that for point vortices the 
variance of the asymptotic Gaussian is = \l N In A^/27r. Figure 1 shows how PNiu/a]^) 
asymptotically approaches a normal distribution. Due to symmetry, Pn{u) = Pn{v) and 
these two distributions have been combined into a single pdf. Further, pat is even and we 
combine positive and negative velocities. A least-squares fit of the small velocity portion 
of the distribution to a Gaussian provides a variance of approximately one (as expected 
for —>■ oo) and an amplitude 0.9 of that of a normal distribution. One sees that even 
though the central limit theorem formally applies, the convergence is extremely slow and 
even for large N, has significant high- velocity tails. Further, the small velocity part of 
Pat is well approximated by a Gaussian with the asymptotic variance, but with amplitude 
less than one due to the significant fraction of events contained in the high-velocity tails. 
In what follows, we show that both the Gaussian portion of the pdf at small velocities and 
the high-velocity tails play an important role in the advection process. 

In the remainder of this paper we focus on systems of 100 point vortices. Similar 
results are obtained with other values of between about 10 and a thousand. From the 
perspective of dynamical systems, an ensemble of 100 point vortices has a 200-dimensional 
phase space; it is thus quite high dimensional. However, it is still reasonable to numerically 
integrate the equations for significant times. From the perspective of the central limit 
theorem, due to the slow non-normal convergence, = 100 is not large enough for pioo to 
be completely Gaussian. 

We next consider several different initial configurations of the A^ vortices. The 
circulations of the vortex populations are characterized by T, the average of the absolute 
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value of the individual circulations. We focus on configurations with F = 1, and with 
random individual circulations uniformly distributed in the range O.Sr < |rj| < 1.2r. The 
variation in the individual Fj's is included to break any symmetry in vortex pair interactions 
that may arise in a system with identical |rj|'s. For each configuration, the initial vortex 
positions are randomly chosen in the square domain [0, 27r]^. For simplicity we only consider 
initial conditions with = Py = 0. 

The set of initial configurations discussed above can have a range of energies, 
E = if ({xj}) /47r^. Close opposite-sign vortex pairs give a large negative contribution to 
the energy, while close same-sign pairs provide a positive contribution. Thus, depending 
on the random initial positions of the vortices, the energy of the system may take different 
values. In Figure 2 we show the density of states, i.e., the number of states with energy in 
a given interval, obtained from 10^ different realizations of the initial vortex configurations 
described above with N — 100. This distribution has a mean of -0.076, a variance of 1.18, 
and a skewness coefficient of 0.85. In the following we integrate the motion of 300 random 
configurations extracted from this ensemble for a relatively short time, T = 0.1. The two 
vertical dotted lines indicate two arbitrarily chosen configurations that have been integrated 
for a much longer time, T = 30. 

The relationship between vortex speed and nearest neighbor distance for all the 100 
vortices in the 300 randomly chosen configurations is seen in Figure 3. Large speeds are 
associated with close nearest neighbors. From this we may conclude that the high velocities 
in the tail of the pdfs are due to the presence of a single close vortex and not to the 
superposition of the contributions of several different vortices. Thus, high-velocity tails in 
Pioo are essentially two-vortex phenomena and not many-body effects. This implies that in 
the limit — > cxd, as the tails of pi\r disappear, the importance of two- vortex interactions 
goes to zero. However, given the slow convergence of pn, two-vortex interactions will be 
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important even for very large N. 

4. Long-Time Integrations 

In this section we focus on the phase-space trajectories obtained by a relatively 
long-time integration of the two randomly selected initial configurations indicated in 
Figure 2. These initial configurations have N — 100 vortices and Np — 100 passives. The 
configurations were integrated until T = 30 using a fourth-order Runge-Kutta method with 
fixed time step At = 10"^. The positions and velocities of the vortices were saved every 
Ats = 10-2. 

The Cartesian velocity pdfs, averaged over time and all vortices, denoted by p{u), 
are obtained by measuring u and v for all vortices every Ats throughout the integration; 
again, from symmetry p{u) — p{v) — p{—u) — p{—v) so ±it and ±v are combined into a 
single pdf. The pdfs for the vortex and passive velocities in each of the two solutions are 
shown in Figure 4. The central part of each pdf, \u\ ^ 20, is well approximated by the 
same Gaussian which fits the central part of pioo! it has the theoretical asymptotic width 
(Tioo = y'lOO In 100/27r ~ 8.56 and amplitude 0.9 of a Gaussian with unit normalization. 
Thus the vortices and passives in the two solutions all have the same small- velocity pdf as 
each other and as an ensemble of random initial conditions. 

The high- velocity tails of p for the two solutions are, however, significantly different. 
One solution has a significant excess over pioo at large velocities, while the other has a 
deficit. Note that both solutions, however, have a deficit at extremely large velocities, 
u ^ 250 ~ 30(7100 • The pdf with smaller tail in the vortices comes from the same integration 
as the pdf with the larger tail in the passives, and vice-versa. This is not significant, and is 
merely the result of picking only two random initial conditions. 
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The lack of convergence of the pdfs over the time of the integration T indicates that 
T is not long enough for the system to forget its initial conditions. By other measures, 
however, T is very long. A typical eddy turnover time, estimated as the time for a vortex 
at a typical vortex separation 27r/v^ from another vortex to rotate in a complete circle, 

1 /2 

is approximately 0.4. Another measure of the eddy turnover time is (Z^Tf) ~ 0.1. Thus 
T is O(IOO) eddy turnover times. A completely different measure of the length of the 
integration is that in a time of T the particles travel around the basic periodic domain 
several times. All of this indicates that the correlation time for the tail of the pdf is very 
long, and that the point- vortex system displays a very long memory with an associated lack 
of ergodicity over large, if not infinite, time intervals. While this is consistent with what 
has already been observed for systems of a few point vortices, it is somewhat surprising 
that this phenomenon persists in a system with so many degrees of freedom. 

A standard statistic for studying velocity time series is the Lagrangian velocity 
autocorrelation -R(r), 

_ (u(t)-u(t + r)) 

where ( ■ ■ ■ ) represents an average over time and over particles. The autocorrelations from 
the long phase-space trajectories for the two populations of vortices and passives are shown 
in Figure 5. As in the pdfs, there is no significant difference between the vortices and 
passives. After a very short period with steep decay, the autocorrelation of the vortices 
and passives follows an exponential extremely well. For r greater than those shown in the 
Figure, R{t) oscillates around zero. The Lagrangian decorrelation time of the phase-space 
trajectories, T^, is estimated by fitting Rij) ~ exp(— r/T^) for the period r = 0.02 to 
r = 0.2, giving Tl = 0.09 ± 0.01. A theoretical estimate of the order of magnitude of Tl 
can easily be obtained by assuming that is the time it takes for a vortex moving at a 



typical velocity (Jn = \/N\nN/2-iT to cross a distance equal to the typical vortex separation 



2-k/^/N; this estimate gives Tl ~ (2tt)^^'^/N\/\iiN ~ 0.07, quite close to that effectively 
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observed. Note also that the common method of estimating the Lagrangian integral time 
scale, Tj = i?(r)rfr, provides misleading results because the very long memory of the 
system implies that extremely long integration times are needed to reliably calculate Tj. 

Individual vortex and passive trajectories, some examples of which are shown in 
Figure 6, are in general quite complex. The first two panels of Figure 6 show examples of 
an intermittently fast moving vortex and of a vortex randomly selected from the population 
of vortices with moderate average speed. The single fast vortex shown in Figure 6(a) 
comprises 33% of the extreme tail (velocity greater than lOaioo) of the velocity pdf p for 
its integration. Figure 6(c) shows the particle trajectory of a randomly selected passive 
particle. The particle trajectories of the slow vortex and the passive are qualitatively 
similar, although the passive does have more tight loops, indicating that at those times it is 
close to a vortex. The particle trajectory of the fast vortex is qualitatively different than 
the other two. 

The time series of the speed of the fast vortex. Figure 7, shows clear bimodal behavior 
with the vortex jumping between slow and fast episodes. A time series of the distance to 
the closest vortex and its identity. Figure 8, shows that episodes of extremely fast motion 
coincide with the vortex being close to another, oppositely-signed, vortex, i.e., it is a 
member of a close dipole. These episodes start and end with changes in the identity of the 
close vortex. Furthermore, the three separate fast episodes are separated by extremely short 
periods of slow motion where the closest vortex is significantly further away and changes 
identity rapidly. The ratio of the distance of the second closest vortex to the distance of 
the closest vortex. Figure 9, demonstrates that the fast episodes are primarily the result 
of close dipoles. However, note that during the period from t ^ 22 to t ^ 24: this distance 
ratio is smaller than during other fast periods (but still larger than during slow periods), 
and there are small oscillations in the nearest neighbor distance (Figure 8). Analysis of the 
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individual vortex trajectories shows that during this period the second closest vortex has 
the same identity and there is thus a temporary three-vortex bound state. 

Before closing this section, we note that we have not detected any significant difference 
between vortex and passive statistics, in the sense that the internal variability in each of 
these populations is by far larger than any difference between vortices and passives. This 
is not too surprising if we consider that the contribution of each individual vortex to the 
total vorticity field is rather small for = 100, and each individual vortex in the collective 
field of the other vortices behaves similar to a passive. Significant differences do exist, 
however, in strong two-body interactions: close vortices accelerate each other and may 
undergo episodes of fast displacement, while passives trapped in the vicinity of a vortex do 
not infiuence its motion. 

When passives are close to a vortex, they spin with large velocity and create 
high-velocity components in the pdf. Moreover, passives close to a vortex tend to remain 
associated with it for very long times, both on the infinite plane0 and in periodic domains.0 
The presence or absence of high- velocity "bumps" in the pdf is thus entirely determined by 
the initial conditions, i.e., whether or not the system happens to start with close passives. 
Only by integrating for times much longer than the trapping time is it possible for the 
time average pdf to converge to the ensemble average pdf, an issue which is related to the 
possible non-ergodicity of passive particle motion in point vortex systems.0 At present, it 
is unclear whether the trapping islands around the vortices have a time-asymptotic nature 
or disappear at finite (but long) time. Preliminary runs of few- vortex systems^ have shown 
that the trapping islands exist up to at least t = 3000. By comparison, bound states of 
vortex couples tend to live for much shorter times. 

When a passive is close to a vortex, it thus behaves similar to the vortex with which 
it is linked, with an additional fast rotational component. Whenever the vortex undergoes 
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a fast displacement due to vortex dipolar coupling, the trapped passive does as well. A 
detailed theoretical comparison of the long-time displacement statistics of vortices and 
passives should thus carefully evaluate the probability that a passive be trapped near a 
vortex (which is related to the initial conditions and determined mainly by geometric 
factors) and the lifetime of vortex-vortex and vortex-passive couples. On a purely heuristic 
basis, the present simulations have shown no detectable differences between the gross 
statistical behavior of vortices and passives, at least on the time scales we have considered. 



5. Properties of single-particle dispersion 

In this section we study the properties of single-particle dispersion in point-vortex 
systems. A classic stochastic model for describing the motion of advected particles is the 
Ornstein-Uhlenbeck (OU) processB 

dx = udt, du = ——dt H ^(iW, (6) 

where W is the Wiener process, and dW is a Gaussian random increment with (dW) = 
and {dWait) dWjjit')) = 2 5a^p5{t — t') dt, where ( ■ • ■ ) indicates an average over an 
ensemble of independent realizations, and Greek subscripts indicate vector components. 
The OU process is characterized by two parameters: a velocity scale (Tqu which determines 
the variance of the velocity pdf, and a timescale which is the Lagrangian decorrelation 
time of the exponential velocity autocorrelation Rouij) = exp(— r/T^). Recently, the OU 
process has been used to model single-particle dispersion in the atmosphere and the ocean.!! 

Here, the OU process is a very natural choice for a stochastic model. As we have 
already seen, the velocity pdf has a significant Gaussian component, leading us to pick a 
Gaussian process. Additionally, the Lagrangian velocity autocorrelation has a period of 
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exponential decay, which leads to the OU process. Furthermore, the two parameters of the 
OU model have already been determined: Oou = c"ioo ~ 8-56 and Tl = 0.09. The major 
failure of the OU process is that it does not capture the high-velocity tails of the vortex 
pdfs. A trajectory of an OU particle with these parameters is shown in Figure 6(d). The 
large scale motion is similar to that of a slow vortex or passive, but the small scale motion 
has more sharp turns instead of small loops. 

The OU process may thus give a framework for analyzing the time evolution of the 
vortex statistics, at least for the slow portion of the vortex population. This is in contrast to 
the application of the central limit theorem in Section 3, which only refers to single instants 
in time. To explore similarities and differences with the OU process, relatively short-time 
integrations, until T = 0.1, were performed on 300 initial conditions chosen randomly from 
those shown in Figure 2. Each initial condition has 100 vortices and 100 passives. This 
dataset is thus the same size as a single long integration. The velocity pdfs averaged over all 
initial conditions, as well as the time and initial condition averaged pdfs, are equal to pioo- 

We now compare the vortex displacements, Ax{t) = x{t) — a;(0) with the displacements 
of the OU process. As is usually done in the study of systems with periodic boundary 
conditions, we compute displacements from the unfolded trajectories, i.e. by taking into 
account the number of times a particle has wound around the periodic domain.lli Again, due 
to symmetry, we can combine x and y displacements. The OU process has displacements 
with a Gaussian pdf whose width a^x is 



The pdf of vortex displacements at time delay r can be compared to the OU pdf by scaling 
each vortex displacement by the above o"ax(t)- The resulting set of pdfs, as obtained by 
considering all 100 vortices for all 300 initial configurations, are shown in Figure 10 for 
time delays r = 0.01,0.02, ...,0.09. This figure shows that for small displacements the pdfs 




(7) 



-16- 



agree with the OU result. However, there are significant large-displacement tails, which are 
largest for small times, and decrease as r increases. 

To identify the origin of the large-displacement tails, in Figure 11 we show the 
displacements at r = 0.05 versus the nearest neighbor distance, for both same-sign and 
opposite-sign neighbors. The large displacements are associated with close opposite-sign 
vortices, i.e. a vortex dipole, while there is no correlation between large displacements and 
close same-sign vortices. Further, the value of the large displacement is well predicted by 
the motion of an isolated dipole. 

We can conclude that the small displacements are essentially a mean-field process, in 
that 1) they are due to the action of many other vortices and 2) they are well captured 
by a stochastic model such as the OU process, which is known to work well for extremely 
high-dimensional systems such as Brownian motion. On the other hand, the large 
displacements are essentially due to two-body interactions, namely, temporary coupling of 
opposite-sign vortices. Same-sign couples rotate rapidly about their center of vorticity, and 
otherwise behave similarly to a single vortex with larger circulation. Thus, they contribute 
to the high-velocity tails, but do not contribute to the large-displacement tails. Similar 
results are observed for the passive particles, with the usual caveat that when passives 
couple with a vortex they do not affect its velocity, but they can be trapped in the vicinity 
of a vortex undergoing fast dipolar motion. 

We next turn to the temporal evolution of the average vortex and passive displacement, 
i.e, to absolute or single-particle dispersion. The single-particle dispersion of a system of 
moving particles is defined as 

A\T,t)^{{^{t + T)-^{t)f)^, (8) 

where { ■ ■ ■ )p indicates an average over the ensemble of particles. When the particle motion 
is statistically stationary, as is the case here, becomes independent of t. The finite-time 
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dispersion coefficient is defined as -D(t) = A^/2r. For motions in differentiable velocity 
fields, the short time dispersion is ballistic, i.e. A'^{t) oc for t << Tj where Tj is the 
Lagrangian integral time scale, defined as the integral of the velocity autocorrelation. If 
Tf is finite, then at large times the dispersion becomes Brownian-like, i.e. A'^{t) ~ r. The 
usual dispersion coefficient is thus defined in the infinite time limit, D = lim^^oo D[t). 

For the OU process, the appropriate ensemble of particles is a collection of independent 
realizations of the Wiener process. The diffusion coefficient then becomes 



which has the infinite time limit Dou = 2o"o„Ti. 

For the point-vortex system, the appropriate ensemble in Equation (^) is the collection 
of vortices or passives at a single time in a single integration. As for displacements, 
dispersion is calculated using unfolded trajectories. The four resulting diffusion coefficients 
are shown in Figure 12. There is no significant difference between the vortices and passives, 
and thus the best estimate is obtained by averaging the four diffusion coefficients. The 
diffusion coefficient for the OU process. Equation (^, is shown for comparison. Also shown 
is the expected error of the mean for four ensembles of 100 OU particles, which is calculated 
with standard statistical methods. For short times, t ^ 0.1 ~ T^, the OU process captures 
the dispersion reasonably well. For intermediate times, <^t <^2 ^OTl, the OU process 
overestimates the dispersion. For long times, t ^ 20Tl, the diffusion coefficient appears to 
grow beyond that of the OU prediction. This anomalously large single-particle dispersion is 
presumably due to the tails in the velocity pdfs and large displacements of long-lived close 
dipoles. It is an interesting question, which remains unanswered at the moment, whether 
the single-particle dispersion becomes Brownian at very large times in point- vortex systems. 
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6. Conclusions 

In this paper we analyze the Lagrangian dynamics of vortices and passives in 
high-dimensional, N = O(IOO), point- vortex systems. We find that the behavior can 
be understood by partitioning the system into two components: (1) mean-field behavior 
resulting from the collection of distant vortices, and (2) high-velocity behavior resulting 
from close pairs (and, occasionally, triplets or more). 

The partition in the statistical dynamics is due to the slow, non-normal convergence 
of pn- In the limit N oo the velocity pdf is Gaussian and only mean-field behavior 
occurs. The mean-field behavior is well modeled as an Ornstein-Uhlenbeck (OU) stochastic 
process. This part of the dynamics has a short correlation time, = 0.09, which can be 
estimated by simple arguments based on the turn-over times for average vortex separation. 
However, due to the l/|u|^ behavior of the velocity pdf for a single vortex at large |u|, the 
convergence to the Gaussian implied by the central limit theorem is extremely slow with 
N, and the variance of the Gaussian component scales as y/N\nN rather than the more 
common 

The slow convergence is important for physically relevant, large but finite, values of A^, 
where there is a significant, non-Gaussian, high- velocity tail in the pdf. This high- velocity 
tail has a very long correlation time. Close pairs typically last for much longer than T^; 
indeed, they can last so long that even our longest integrations, of O(300Ti), do not 
yield reliable statistics for the distribution of pair lifetimes. These long-lived pairs exhibit 
nonergodic behavior for at least the time of our simulations, which is consistent with 
previous evidence for long-time nonergodicity with only a few point vortices.0' 

As a result of the high-velocity tails, the OU process does not accurately capture the 
single-particle dispersion and the pdf of particle displacements. The small displacements 
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scale approximately as in the OU model, but the large displacements show no such scaling 
due to close dipoles. At very short times, of course, the displacements exhibit ballistic 
scaling, with D{t) oc r. After the first significant loss of velocity correlation, at t ~ 0{Tl), 
the OU model captures the average dispersion reasonably well. At later times, t ^ 1 - 5, 
the OU model overestimates the single-particle dispersion, and at very large times, t > 20, 
the OU model appears to underestimate the dispersion. This suggests that the diffusion 
coefficient may continue to increase resulting in long-time anomalous dispersion, rather 
than saturate as in the OU model around a time of O{10Tl). 

It may be possible to construct a new stochastic process that includes both the 
Gaussian slow vortices and the long-lived, close, high-velocity pairs. The slow vortices 
would be modeled as an OU process. Since the close same-sign pairs rotate in place, these 
vortices should be modeled by a process that has negative velocity correlations at very short 
times. The close dipoles travel large distances, and could be modeled by a Levy walk.0 
In previous studies of point vortices, Viecelli found anomalous relative dispersion at short 
times, which he explained as Levy walks. 

In order to create a combined stochastic process one would need information about the 
close-pair dynamics. This information could be condensed into two joint pdfs and two scalar 
probabilities: the joint pdfs of close-pair velocities and lifetimes, separately for the dipoles 
and same-sign pairs, and the transition probabilities for switching from OU motion to close 
dipoles and same-sign pairs, i.e., the probabilities for close-pair formation. We already 
know that the pdf for close pair velocities, which is related to the joint pdf by integrating 
over lifetime, scales as u^^ for large u. If we are only concerned with displacements, then 
the same-sign pairs can probably be absorbed into the OU component, since their center of 
vorticity moves as a single vortex. If we wish to model the velocities, however, then the fast 
rotation of the same-sign pairs is important. 
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With this information one could construct a stochastic algorithm in which particles 
move with OU random walks, and, at particular instants, have some probability of becoming 
members of a dipole or same-sign pair. Once a member of a dipole, the particle moves for 
some time at a constant velocity in a random direction, with the time and velocity chosen 
from the joint pdf. Once a member of a same-sign pair, the particle has a high oscillating 
velocity for some time, plus a slower OU random-walk velocity, again with the time and 
velocity chosen from the joint pdf. Because of the long correlation time of the close pairs, 
obtaining a good estimate of these joint pdfs would require significantly longer integrations 
than we have done here; therefore, we will not propose here an explicit form for the new 
stochastic model. 

How relevant is the range of N we study here for geophysical flows? The number 
of coherent vortices is not well known in nature, since they are often hard to detect by 
conventional measurement techniques; for example, the vorticity or potential vorticity fields 
are usually poorly sampled. However, we can estimate an upper bound on the number of 
coherent vortices in the large-scale geostrophic regimes of the atmosphere and ocean as the 
ratio of total area to a typical vortex area, neglecting any stacking in the vertical. (We 
do not consider the ageostrophic, smaller-scale regimes, for which 2D vortex dynamics are 
probably less relevant.) A typical geostrophic vortex size is given by the Rossby radius of 
deformation, which, at midlatitudes, is approximately 30 km in the ocean, and 1000 km in 
the atmosphere.S The total areas of the ocean and atmosphere are about 3.5 x 10^ km^ and 
5 X 10^ km^, respectively. Thus, if the area were filled with closely packed vortices, which 
is certainly not the case, the ocean would have O(IO^) vortices and the atmosphere would 
have only O(IO^) vortices. From Figure 1 we see that with this upper bound even the ocean 
would not be in the infinite N regime and the tails of the pdf would be important. Thus, 
the behavior described in this paper is representative of what we expect from the number 
of coherent vortices in the atmosphere and ocean. 
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The point-vortex system is only an approximation to the dynamics of coherent vortices 
in geophysical flows. True coherent vortices have a finite core size and a finite lifetime, 
even though the former is often small and the latter large compared to many geophysical 
phenomena. The existence of a finite core does not significantly affect the dynamics of 
passive tracers, as shown by simulations of two-dimensional turbulence, omce 
extended coherent vortices effectively retain the matter in and near their cores and can 
travel long distances over their long lifetimes, they can have considerable influence over 
large-scale material transports, similarly to what happens for point vortices. On the 
other hand, the dynamics of the vortices themselves may be affected by the presence of a 
finite core. In particular, same-sign vortices can merge, a process which is not included 
in point-vortex dynamics. Although a full answer to the impact of finite cores can only 
come from simulations of extended vortices, we note that many of the specific Lagrangian 
properties discussed here are mainly determined by the behavior of close dipoles. These 
still exist, with similar properties, in the case of extended vortices, suggesting that the 
results found here for point vortices which are due to dipoles may be more general. Future 
studies will address the properties of Lagrangian transport in systems of coherent vortices 
with extended cores and in punctuated point vortex dynamics where instantaneous merger 
is allowed. 

In this paper we have demonstrated several differences of vortex systems from random 
walking. On the other hand, aspects of the overall dispersion behavior seen in Figure 12 
are grossly captured by simple OU diffusion on times very large compared to T^. Thus, 
we see this study as giving some further degree of support for the common practice of 
parameterizing large-scale transports as eddy diffusion. In this context, the eddy diffusion 
coefficient is a function of the coherent vortex population obtained from Equation (P) with 
Gou replaced by o^^ D = NTlItliN/tt. Thus, one may be able to obtain a time-varying 
eddy diffusion by observing how the vortex population changes over time. However, it may 
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eventually turn out that the quantitatively significant differences between vortex dispersion 
and stochastic diffusion imply a significant qualitative difference as well. For example, if 
the diffusion coefficient continues to grow with time, then for times much longer than these 
simulations, eddy diffusion will fail. 
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A. Appendix: The Central Limit Theorem for Point Vortices 

In this appendix we show that in the limit N oo the velocity probability density 
function (pdf) for the velocity due to N randomly placed vortices, p^, becomes Gaussian 



with variance = yAln A/27r. The calculation is based on theorems regarding random 
variables.0 The fact that becomes Gaussian was mentioned previously by Min, et al;@ 
the value of the variance is, to our knowledge, new. 

The relevant theorems discuss the behavior of zn, the sum of N independent identically 
distributed random variables Xj with pdf F{x), scaled by a factor a^: 

Xi + . . . + Xn .... 

zn = • (Al) 
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If in the limit N oo the pdf for zn converges weakly to some distribution G, then F 
is said to be in the domain of attraction of G. The family of distributions which have 
non-empty domains of attraction are called Levy distributions and are identified by an 
index a,0<a<2;a = 2 corresponds to the normal distribution, and a = 1 is the Cauchy 
distribution.!!' S Here we shall only be concerned with a = 2. 



Ibragimov and LinnikB^ prove that F belongs to the domain of attraction of the normal 
distribution if and only if the truncated second moment 



J(X) = x^Fix) dx (A2) 
J-x 

is slowly varying in the sense that 

J"!Lf§-^ (A3) 

for all t. They further show that the scaling of the sum, ctat, is given by the requirement 
that 

hm \ — '- = 1, A4 



for some e > 0. As written. Equation (|A4|) assumes that F has zero mean which is the case 



of interest here; if F has nonzero mean then another term is needed in the equation. 

Here, pi{u) takes the role of F{x), where pi{u) is the pdf for the x-component of the 
velocity, u, produced at the origin from a point vortex with unit circulation placed at a 
random position in the periodic domain [— 7r,7r]^. Since the velocity is a vector u = {u,v) 
described by a joint distribution p{u,v), pi is obtained from 

/oo 
p{u,v)dv. (A5) 
-oo 

Due to symmetry the pdfs for u and v are identical and we need only consider one of them. 
It is sufficient to only consider vortices with positive circulation since changing the sign of 
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the circulation is equivalent to changing the sign of x. Thus, allowing random positive and 
negative unit circulations would not affect pi or p^- 

We are interested in the velocity component u due to vortices, with pdf pn{u), 
where u is the unsealed sum of the velocity components of the individual vortices. Thus, if 
the sum scaled by has a normal distribution, then pjv is Gaussian with variance a^. 

While the joint distribution p{u, v) is rather difficult to determine in detail, changing 
variables allows us to easily integrate over the distribution. The random velocity u is due to 
a vortex at a random position x, so p{u, v) du dv = p{x, y) dx dy. Further, since the vortex 
is placed randomly from a uniform distribution, p{x, y) = I/Att"^. A final change of variables 
to polar coordinates (r, 6) gives 



To determine the bounds of the integral we need the curve r{6,X) such that 
u{r, 6) = X, with v arbitrary. Since we are interested in the behavior of I{X) in the limit 



original integral is over \u\ < X, in polar coordinates the integration is over r > 2sin6'/X; 
this corresponds to integrating over all points outside two circles with radius 1/X centered 
at ±1/X. Along the x-axis, 6* = or tt, r is allowed to go to zero; thus, v = 2cos6/r can 
reach both ±oo and r > 2sin6'/X does not restrict v. Defining R{6) to be the distance 
from the origin to the edge of the periodic domain results in 



Despite the fact that the integral covers small u, the limiting behavior of / only 




(A6) 



of large X we can use the asymptotic form of Equation (|), u{r, 0) 



2sin^/r. Since the 




(A7) 
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depends on the asymptotic form of u{r, 9). Writing the velocity as 



u{r,e) = '^^ + g{r,e), (A8) 
r 

where, from Equation @, (7 ~ 0(r) for small r, gives 

/(X) = -lnX + 0(l). (A9) 

TT 

It is now straightforward to verify that /(X) is slowly varying in the sense of (|X3|), and 
thus pi is in the domain of attraction of the normal distribution. Further, the requirement 
( |A^ ) is satisfied by 

= (AlO) 

Thus, in the limit X ^ 00, is Gaussian with variance a^- Because this result only 
depends on the explicit form of u for small separations, it is independent of boundary 
conditions. 
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Figure Captions 

Figure 1. Probability density function p^iu/uj^) for the instantaneous velocity u 
produced by N randomly placed vortices, scaled by = \J N \\iN/2ti. Solid lines are pn 
for = 10^, 10^, and 10^. Dashed line is a Gaussian pdf with unit variance and amplitude 
0.9. The inset shows same curves on log-log axes, along with the dotted line ~ 1/v?. 

Figure 2. Density of states as a function of energy for 10^ random configurations of 100 
vortices. Dashed lines indicate two arbitrary initial condition used for long integrations. 

Figure 3. Scatter-plot of instantaneous vortex speed versus nearest neighbor distance 
for all vortices in 300 randomly chosen configurations. The dashed lines indicate the speed 
induced by a single vortex with the smallest and largest circulations in the population. 

Figure 4. The velocity pdf averaged over time and particles, p{u), for a) the vortices, 
and b) the passives. In each panel, the two solid lines are p{u) for the two long integrations 
discussed in the text, the dotted line is a Gaussian with variance aioo and amphtude 0.9 
and the dashed line is pioo(w), as in Figure 1, but with the same sample size as p. 

Figure 5. Velocity autocorrelation function R{t), Equation (^), for the two long 
integrations. Solid lines are the vortices, dashed lines are the passives. The dotted line is 
the exponential which best fits the average of the vortices and passives, and has a decay 
time of Tl = 0.09. 

Figure 6. Individual particle trajectories from the long integrations: a) a fast vortex, 
where the period of fast motion, |u| > 70, is shown in bold, b) a slow vortex, c) a slow 
passive, and d) a stochastic particle governed by the Ornstein-Uhlenbeck process. Equation 
(^). The trajectory of the fast vortex appears not to be smooth because of the finite 
plotting interval, Atg. 
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Figure 7. Time series for the speed of the fast vortex shown in Figure 6(a). 

Figure 8. Nearest neighbor distance for the fast vortex shown in Figure 6(a) (lower 
curve). The identity of the nearest neighbor is shown in the upper curve: a hne at 1.3 
indicates a same-sign nearest neighbor, a hne at 1.1 indicates an opposite-sign nearest 
neighbor, and the hne jumps to 1.2 every time the identity of the nearest neighbor changes. 

Figure 9. Ratio of the distances from the fast vortex shown in Figure 6(a) to its second 
and first nearest neighbors. 

Figure 10. Probabihty density function of displacements Ax(r) scaled by the OU 
expected displacement (Tax, Equation (0), for r = 0.1, 0.2, . . . , 0.9, from the 300 short 
integrations. The dashed line is the OU Gaussian. The inset shows the tails for r = 0.1 
(upper solid curve) and 0.9 (lower solid curve). 

Figure 11. Scatter-plot of displacement Ax{t), t = 0.05, versus average nearest 
neighbor distance over the time period for a) opposite-sign pairs, and b) same-sign pairs in 
the 300 short integrations. The solid line in a) is the displacement of a single dipole on the 
infinite plane with |r| = 1. 

Figure 12. Diffusion coefficient -D(r) from the two long integrations for the vortices 
(thin solid curves), passives (thin dashed curves), and the average of the four (thick solid 
curve). Also shown is the OU diffusion coefficient (thick solid smooth curve) and the 
envelope of expected error for four ensembles of 100 OU particles (thick dashed curves). 
The inset shows the same curves for short times. 
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